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FLUCTUATIONS FOR ANALYTIC TEST FUNCTIONS IN THE 

SINGLE RING THEOREM 

FLORENT BENAYCH-GEORGES AND JEAN ROGHET 


Abstract. We consider a non-Hermitian random matrix A whose distribution is in¬ 
variant under the left and right actions of the unitary group. The so-called Single Ring 
Theorem, proved by Guionnet, Krishnapur and Zeitouni [31], states that the empirical 
eigenvalue distribution of A converges to a limit measure supported by an annulus S. In 
this text, we establish the convergence in distribution of random variables of the type 
Tr(/(A)M) where / is analytic on S and the Frobenius norm of M has order ^/N. As 
corollaries, we obtain central limit theorems for linear spectral statistics of A (for ana¬ 
lytic test functions) and for finite rank projections of /(A) (like matrix entries). As an 
application, we locate outliers in multiplicative perturbations of A. 


1. Introduction 

The Single Ring Theorem, by Guionnet, Krishnapur and Zeitouni [31], describes the em¬ 
pirical distribution of the eigenvalues of a large generic matrix with prescribed singular 
values, i.e. an. N x N matrix of the form A = UTV, with U, V some independent Haar- 
distributed unitary matrices and T a deterministic matrix whose singular values are the 
ones prescribed. More precisely, under some technical hypotheses, as the dimension N 
tends to inhnity, if the empirical distribution of the singular values of A converges to a 
compactly supported limit measure 0 on the real line, then the empirical eigenvalues dis¬ 
tribution of A converges to a limit measure /i on the complex plane which depends only 
on 0. The limit measure /i (see Figure 1(a)) is rotationally invariant in C and its support 
is the annulus S' := { 2 ; G C; a < |^| < b}, with a, 6 > 0 such that 



In this text, we consider such a matrix A and we study (Theorem 2.4) the joint weak 
convergence, as N ^ 00 , of random variables of the type 

Tr(/(A)M), 
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for / an analytic fnnction on the annnlns S whose Lament series expansion has nnll con¬ 
stant term and M a deterministic N x N matrix satisfying some limit conditions. These 
limit conditions (see (2)) allow to consider both: 

- flnctnations, aronnd their limits as predicted by the Single Ring Theorem, of linear 
spectral statistics of A (for M = I): 

N 

Tr/(A) = ^/(Ai), 

i=l 

where Ai,..., Aat denote the eigenvalnes of A, 

- hnite ra nk projections of /(A) (for M = y/N x (a matrix with bonnded rank)), like 
the matrix entries of /(A). 

Let ns present both of these directions with more details. 

1.1. Linear spectral statistics of A. As far as Hermitian random matrices are con¬ 
cerned, linear spectral statistics flnctnations nsnally come right after the macroscopic be¬ 
havior, with the microscopic one, in the natnral qnestions that arise (see e.g., among the 
wide literatnre on the snbject, [34, 37, 54, 33, 6, 5, 39, 53, 1, 39, 4, 22, 13]). For nni- 
tary or orthogonal matrices, also, many results have been proved (see e.g. the results of 
Diaconis et al in [27, 28], the ones of Soshnikov in [55] or the ones of Levy and Mai'da 
in [38]). For non-Hermitian matrices, established results are way less numerous: the hrst 
one was [51], by Rider and Silverstein, for analytic test functions of matrices with i.i.d. 
entries, then came the paper [52] by Rider and Virag, who managed, thanks to the explicit 
determinantal structure of the correlation functions of the Ginibre ensemble, to study the 
fluctuations of linear spectral statistics of such matrices for test functions. Recently, in 
[46], O’Rourke and Renfrew studied the fluctuations of linear spectral statistics of elliptic 
matrices for analytic test functions and, in [36, 21], Cebron and Kemp used a dynamical 
approach to study such fluctations on Gh^. The reason why, except for the breakthrough 
of Rider and Virag in [52], many results are limited to analytic test functions is that when 
non-normal matrices are concerned, functional calculus makes sense only for analytic func¬ 
tions: if one denotes by Ai,..., Aat the eigenvalues of a non-Hermitian matrix A, one can 
estimate numbers Tr A^ or Tr((z — A)~^) only when / is analytic. 

For a test function /, one relies on the explicit joint distribution of the Aj’s or on Girko’s 
so-called Hermitization technique, which expresses the empirical spectral measure of A as 
the Laplacian of the function x i —> log | det( 2 ; — A)| (see e.g. [30, 18]). This is a way more 
difficult problem, which we consider in a forthcoming project. 

In this text, as a corollary of our main theorem, we prove that for A = UTV a.n N x N 
matrix of the type introduced above and / an analytic function on a neighborhood of the 
limit support S of the empirical eigenvalue distribution of A, the random variable 

Tr/(A)-E/(A) 
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converges in distribntion, as iV —)■ cx), to a centered complex Ganssian random variable 
with a given covariance matrix (see Corollary 2.8). This is a hrst step in the stndy of 
the noise in the Single Ring Theorem. We notice a qnite common fact in random matrix 
theory: the random variable 

N 

Tr/(A) -ETr/(A) = ^/(A.) - E/(A,) 

i=l 

does not need to be renormalized to have a limit in distribntion, which reflects the eigen- 
valne repulsion phenomenon (indeed, would the Aj’s have been i.i.d., this random variable 
would have had order a//V). 

Next, two corollaries are given (Corollaries 2.12 and 2.14), one about the Bergman kernel 
and the resolvant and one about the log-correlated limit distribution of the characteristic 
polynomial out of the support. 

It should be noted that the class of test functions studied - / analytic on a neighborhood 
of the annulus where the eigenvalues Xi locate asymptotically - is not rich enough to fully 
characterize the fluctuations the spectrum. For example, not all smooth functions on the 
annulus can be approximated by analytic functions. Thus while these results do give insight 
into the fluctuations, the full study of the fluctuations would have to go beyond the realm 
of analytic test functions. 

1.2. Finite rank projections and matrix entries. A century ago, in 1906, Emile Borel 
proved in [19] that, for a uniformly distributed point (Xi,... ,Xn) on the unit Euclidian 
sphere the scaled first coordinate a/XXi converges weakly to the Gaussian distribu¬ 

tion as the dimension N tends to inhnity. As explained in the introduction of the paper 
[3] of Diaconis et ai, this means that the features of the “microcanonical” ensemble in a 
certain model for statistical mechanics (uniform measure on the sphere) are captured by 
the “canonical” ensemble (Gaussian measure). Since then, a long list of further-reaching 
results about the asymptotic normality of entries of random orthogonal or unitary matrices 
have been obtained (see e.g. [3, 40, 23, 26, 35, 10, 12]). 

In this text, as a corollary of our main theorem, we prove that for A = UTV an X x 
N matrix of the type introduced above, f{z) = an analytic function on a 

neighborhood of the limit support S of the empirical eigenvalue distribution of A and a, b 
some unit column vectors, the random variables of the type 

\/X(bV(A)a-aob*a), 

which can be seen as particular cases of the variables of the type Tr(/(A)M)—E Tr(/(A)M) 
for M = V^ab*, converge in joint distribution, as X —>■ oo, to centered complex Gauss¬ 
ian random variables with a given covariance matrix (see Gorollary 2.16). This allows for 
example to consider matrix entries of /(A), in the vein of the works of Soshnikov et al. 
for Wigner matrices in [47, 50] (see Gorollary 2.19 and Remark 2.20). It also applies to 
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the study of finite rank perturbations of A of multiplicative type: the BBP phase transi¬ 
tion (named after the authors of the seminal paper [7]) is well understood for additive or 
multiplicative perturbations (A = A -|- P or A = A(I -I- P)) of general Hermitian models 
(see [48, 20, 14] or [7, 15]), for additive perturbations of various non-Hermitian models 
(see [56, 16, 45, 17]), but multiplicative perturbations of non-Hermitian models were so 
far unexplored. In Remark 2.17 and Figure 1, we explain briefly how our results allow to 
enlighten a BBP transition for such perturbations. 

1.3. Organisation of the paper and proofs. In the next section, we state our main 
theorem (Theorem 2.4) and its corollaries. The rest of the paper is devoted to the proof 
of Theorem 2.4, to the proof of Corollary 2.14 and to the appendix. 

Theorem 2.4 is proved in three steps. First, we do a cut-off approximation to replace the 
analytic functions / in the random variables Tr(/(A)M) by polynomials. The estimation 
of the error term in this cut-off is far from obvious relies on non-asymptotic estimates from 
[11] and [16]. Then, to prove the convergence to Gaussian random variables, we perform a 
moment calculation, using Weingarten calculus (for the asymptotic hne moments of Haar- 
distributed unitary matrices). Weingarten calculus is the theory, due essentially to Collins 
and Sniady, of the joint moments of entries of Haar-distributed unitary matrices. We 
summarize the necessary ingredients of the theory in Appendix A. The third step is the 
computation of the limit covariance. 


2. Main result 

Let A be a random N x N matrix implicitly depending on N such that A = UTV, 
with U, V, T independent and U, V Haar-distributed on the unitary group. We make the 
following hypotheses on T; 

Assumption 2 . 1 . As A —)• oo, the sequence (A“^ Tr TT*)^'^^ converges in probability to 
a deterministic limit b > 0 and there is M < oo such that with probability tending to one, 
l|T||op < M. 

Assumption 2 . 2 . With the convention 1/cx) = 0 and 1/0 = oo, the sequence 

converges in probability to a deterministic limit a > 0. If a > 0, we also suppose that there 
is M' < oo such that with probability tending to one, ||T~^||op < M'. 

The following, seemingly purely technical, assumption, which could possibly be relaxed 
following Basak and Dembo’s approach of [8], is made to control tails of Laurent series but 
can be removed if the //s have hnite Laurent expansion, like in Corollary 2.9 or in Remark 
2.20. Precisely, we need it to cite some estimates from [32], where they were proved under 
this assumption. 
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Assumption 2.3. There exist a constant k > 0 snch that 

Ini(z) > n-^ => A-MlniTr((^-VTT*)-M| < 

For / an analytic fnnction on a neighborhood of the annnlns 

S' := { 2 ; G C; a < | 2 ;| < b}, 

the matrix /(A) is well dehned with probability tending to one as N ^ 00 , as it was proved 
in [32, 11 ] that the spectrnm of A is contained in any neighborhood of S with probability 
tending to one. We denote the Lanrent series expansion, on S, of any snch fnnction / by 

Theorem 2.4. For each N > 1, let Mi,..., be N x N deterministic matrices such 
that for all i,j, as N ^ oc, 

^TrMj —)■ Ti, ^TrMiM* —^ ^TrMiM^ —^ Pij ( 2 ) 

Let /i,..., /fc be analytic on a neighborhood of S. Then, as N ^ 00 , the random vector 

(Tr/,(A)M,-ao(/,)TrM,)' (3) 

\ /j=i 

converges to a centered complex Gaussian vector {G{fi), ■ ■ ■ ,G{fk)) whose distribution is 
defined by 

E G{fi)G{fj) = X] ((r- - ^)riTj + /3ij) {an{fi)a-n{fj) + a-n{fi)an{fj)) 

n>l 

E Q{fi)Q{fj) = X] ((r - + aij) {an{fi)an{fj)b^'^ + a_„(/i)a_„(/j)a“^"). 

n>l 

Remark 2.5. In (3), 

Tr/,(A)M,-a„(/,)TrM, 

rewrites 

Tr/,(A)M,-ETr/,(A)M,. 

Indeed, E/(A) = ool, as a conseqnence of the fact that for any n 7 ^ 0, for any 6 * G M, 
A” e^^A"', which follows from the invariance of the Haar measnre. 

Remark 2.6. Note that if a = 0, as the //s are analytic on S, we have a_„(/j) = 0 for 
each n > 1 and each j, so that the above expression still makes sense. Besides, it seems 
reasonable to verify that the two series above converge: 

V'R|an(/i)||an(/j)|&^” < (max |a„(/j)| 6 ") ^n|an(/i)| 6 " < cx) 

f ^ n>l ^ ^ 

n>l n>l 

E n\an{fi)\\a_n{fj)\ < (max |a„(/i)| 6 ") y^n|a_„(/j)|a“’" < cx). 

n>l ' ^ 

n>l n>l 
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Remark 2.7 (Relation to second order freeness). A theory has been developed recently 
about Gaussian fluctuations (called second order limits) of traces of large random matrices, 
the most emblematic articles in this theory being [41, 42, 43, 24], Theorem 2.4 can be 
compared to some of these results. However, our hypotheses on the matrices we consider 
are of a different nature than the ones of the previously cited papers, since the convergence 
of the non commutative distributions is not required here: our hypotheses are satisfied 
for example by matrices like Mj = a/ZV x(an elementary N x N matrix), which have no 
bounded moments of order higher than two. 

Our two main applications are the case where the Mj’s are all I (Corollaries 2.8 and 2.9) 
and the cases where the M^-’s are y/N times matrices with bounded rank and norm, like 
elementary matrices (Corollaries 2.16 and 2.19). In the case M = I, we immediately obtain 
the following corollary about linear spectral statistics of A. 

Corollary 2.8. Let /i, ... ,/fc be analytic on a neighborhood of S. Then, as N ^ oo, the 
random vector 

(Tr/,(A)-iVao(/,))' 

\ / j=i 

converges to a centered complex Gaussian vector {Q{fi ),..., G{fk)) such that 

mmif,) = E I ^ I n (/j ) 

n>l 


For n > 1, let us define the functions 


Id ■ 

These functions (plus the constant one) define a basis of the space of analytic functions on 
a neighborhood of S and we have the change of basis formula 


a_r 


= ao + ^cfipf{z)+ c^ip^{z) Vn > 1, 

nSZ n>l 

implying that 

EM/)iv”+ia_„(/)iv^” = 2Eiy(/)r+fc(/)i 


,-n h-n 


—a^ ^ ^ ^ ^ ’ 


n>l 


n>l 


Besides, these functions allow to identify the underlying white noise in Theorem 2.4 (we 
only state it here in the case Mj = I, but this of course extends to the case of general 
Mj’s, allowing for example to state analogous results for the matrix entries). 


Corollary 2.9 (Underlying white noise). The finite dimensional marginal distributions of 


(Tr(p+(A))„>iU(Tr(^„(A)) 

n>l 





FLUCTUATIONS FOR ANALYTIC TEST FUNCTIONS IN THE SINGLE RING THEOREM 


7 


converge to the ones of a collection {Qf^)n>i U {Qn )n>l of independent centered complex 
Gaussian random variables satisfying 

EieJl" = 2 i E(ej)" = ±2n(o/6)" 


Remark 2.10 (Ginibre matrices). In the particular case where A is a Ginibre matrix (he. 
with i.i.d. entries with law J\fc{0, N~^)), we reproduce the result of Rider and Silverstein 
[51], noticing that in this case a = 0 and 6 = 1, so that a„(/) = 0 when n < 0 and 
^G{fi)Q{fj) = 0, and, for dm{z) the Lebesgue measure on C, 


1 

71 

1 


'hl<i 


^Mz)-^fjiz)dm(z) 




^ 2|:|<1 4jr 2circle(l+e) 2Circle(l+e) (^1 V (^2 

MMJAS) 1 


J^ircle(l+e) J^ircle(l+£) ^ 


[ nn' 


z\^~^ 

-) [£l dm(.)d{.d& 




' Circle(l+6:) ^Circle(l+e:) 

y'na„(fi)a„{Jj) 


/.Ki)/2(&)5Z’>te^2)'" dfidfe 


n>l 


n>l 


Remark 2.11. If T = I, and the /^’s are polynomial, we reproduce a result of Diaconis 
and Shahshahani [27, Theorem 1] on the limit joint distribution of 

(Tiiu*));*.,, 

where U is Haar-distributed. Actually, the Gorollary 2.8 is slightly stronger, since the 
result holds for A = UT or A = UTV as long as T satisfies 

Mm llV(TT-) = hm 1 T1((TT-)-') = 1, (4) 

N^oo I\ N^oo I\ 

in which case A may be seen as a multiplicative perturbation of U. Indeed, (4) implies 
that all singular values of T are close to 1. The matrix T satishes the condition (4) for 
example if it is diagonal and all its diagonal coefficients are equal to 1 except o{N) of them 
(which stay away from 0 and cx)). 

Corollary 2.12 (Bergman kernel and resolvant). The random process 

(Tr(^-A)“')|,|<, U (Tr(z-A)-i)|^|^^ 

converges, for the finite-dimensional distributions, to a centered complex Gaussian process 


i^z)\z\<a i'^z)\z\>b 


with covariance defined by 

— b^ 
^GzGz' = 


( 6 ^ — zz'y 


EnMz' = 


(a^ — zz'y 


^GzTLz' = 


{z' — zY 















and by the fact that 

ye e M, 
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|2|<a 


u 




n,) 


|2|>6 


law / ^ 

= [yz 


|2|<a 


U in,) 


\z\>b ' 


Remark 2.13. The kernel of the limit Gaussian analytic function, in the previous corollary, 
is, up to a constant factor, the Bergman kernel (see [9, 49]). 


Corollary 2.14 (Characteristic polynomial out of the support). The random process 
(log|det(^- A)| -TrlogT)|^l^^ U (log | det( 2 ; - A)| - Alog |^|)|^l>^ 
converges, for the finite-dimensional distributions, to a centered real Gaussian process 

^^ z )\ z\<a ("^ 2:)|^|>6 

with covariance defined by 


2Eg,g,f = - log 


zz' 



, 2E g,n,^ = - log 


, 2E n,n,i = — log 

1 - ^ 



zz' 



1 - - 


Remark 2.15. As z ^ z' both tend to the same point on the boundary of S', the above 
covariances are equivalent to — log \ z — z'\. In the light of the log-correlation approach to 
the Gaussian Free Field (see [29]), it supports the idea that on the limit support S, the 
characteristic polynomial of A should tend to an object related to the Gaussian Free Field, 
as for Ginibre matrices (see Corollary 2 of [52]). It would be interesting to see to what 
extent such a convergence depends on the hypotheses made on the precise distribution of 
the singular values of T. 


In the case = \/Naijh*, we immediately obtain the following corollary: 

Corollary 2.16. For each N > 1, let ai, bi,..., a^, be deterministic column vectors 

with size N such that for all i,j, as N ^ oo, 

a*aj —^ Kif e C ; b^a^- —^ G C ; h*hj —^ e C (5) 

Let fi,..., fk be analytic on a neighborhood of S. Then, as N ^ oo, the random vector 

Gv(b*/,(A)a, - b*a,.a„(/,))‘_^ (6) 

converges to a centered complex Gaussian vector (^(/i),..., gifk)) such that 
Egifi)gifj) = '^K^fiK)'^{anifi)a-nifj) + a-nifi)anifj)) 

n>l 

^eUi)sUi) = + 

n>l 

Remark 2.17 (Application to multiplicative finite rank perturbations of A). The previous 
corollary has several applications to the study of the outliers of spiked models related to 
the Single Ring Theorem. It allows for example to understand easily, using the techniques 
developed in [16], the impact of multiplicative finite rank perturbations on the spectrum 
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of A (whereas only additive perturbations had been studied so far). For example, one can 
deduce from this corollary that for P a deterministic matrix with bounded operator norm 
and rank one, if one dehnes A := A(I + P) and A := A(I + AP), then 

• the matrix A has no outlier (he. the support of its spectrum still converges to S'), 

• the matrix A has no outlier with modulus > b, but each eigenvalue A of P such 
that |A| > a~^ gives rise to an outlier of A located approximately at —A“^ (besides, 
when the multiplicity of A as an eigenvalue of P is 1, the fluctuations of the outlier 
around — A“^ are Gaussian and with order iV“^/^). 

This phenomenon is illustrated by Figure 1. 



(a) Spectrum of A 



(b) Spectrum of A 

A(I + P) 



(c) Spectrum of A := 
A(I+AP) (small circles are 
centered at the theoretical 
limit locations of the out¬ 
liers) 


Figure 1. Outliers/lack of outliers for multiplicative perturbations: sim¬ 
ulation realized with a single 10^ x 10^ matrix A = UTV when the singu¬ 
lar values of T are uniformly distributed on [0.5,4] and P = diag(—2, (0.8 -|- 
0.5i)“^, 1/3, 0,..., 0). As predicted, none of these matrices has any outlier outside 
the outer circle, nor do the two first ones inside the inner circle, but A has two 
outliers inside the inner circle, close to the predicted locations. 


To state the next corollary, let us hrst give the dehnition of Gaussian elliptic matrices. 

Definition 2.18 (Gaussian elliptic matrices). Let p E C such that |p| < 1. A Gaussian 
elliptic matrix with parameter p is an A x A Gaussian centered complex random matrix 
X = {xij) satisfying : 

(1) the random vectors are independent, 

(2) Vi, E \xii\^ = 1 and Ex^ = p, 

(3) Vi 7 ^ j, E \xij\^ = 1, Ex^j = 0, ¥.XijXji = p and Ex^x/i = 0. 
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This matrix ensemble was introduced by Girko in [30] , and its name is due to the fact that 
its empirical eigenvalue distribution is the uniform distribution inside an ellipse. In the 
case where p = 0 (resp. p = 1), we get a Ginibre (resp. GUE) matrix. 

Corollary 2.19. Let f be analytic on a neighborhood of S such that 

X] |an(/)P&^’" + |a-n(/)Pa"^” = 1- 

n>l 

Let k be a fixed positive integer and let I = I{N) be a (possibly N-dependent) subset of 
{1,..., N} with cardinality k. Let us define the random k x k matrix 

x„ := Gv [/(A)., - a„(/)%]. 

Then, as N ^ oo, the matrix Xat converges in distribution to a k x k Gaussian elliptic 
matrix X with parameter p, 

p := 2^an(/)a_„(/). 
n>l 

Remark 2.20. a) In the case where f{z) = z, we rederive the well-known result that 
any hxed-size principal submatrix of y/NU converges to a Ginibre matrix (see e.g. 
[3, 35]). 

b) By Gorollary 2.19, the statement of the hrst part of this remark happens to stay 
true, up to a constant multiplicative factor, if U is replaced by A = UTV or even 
by A” or by /(A) if / is analytic in a neighborhood of the disc B{0, b). 

c) It also follows from what precedes that for any n > 1, any sequence of principal 
submatrices with hxed size of yjN and yjN— converge in 
distribution to a GUE matrix and i times a GUE matrix, both being independent. 

3. Proof of Theorem 2.4 

To avoid having to treat the cases a > 0 and a = 0 separately all along the proof, we 
shall suppose that a > 0 (the case a = 0 is more simple, as sums run only on n > 0). 
Besides, note that by invariance of the Haar measure, the distribution of the random 
matrix A depends on T only through its singular values, so we shall suppose that T = 
diag(si,..., Sat), with Si > 0. At last, as the limit distributions, in Theorem 2.4, only 
depend on T only through the deterministic parameters a, b, up to a conditioning, one can 
suppose that T is deterministic (and that both ]|T]|op and ||T~^||op are uniformly bounded, 
by Asssumptions 2.1 and 2.2). 

3.1. Randomization of the M^’s. Let us define W := VU. The random matrix W 
is also Haar-distributed and independent from V. Besides, for each j, as A = UTV = 

V*WTV, 

Tr/j(A)M^- = TrV*/i(WT)VMj = Tr/j(WT)VMjV* 
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As a consequence, we shall suppose that A = UT (instead of A = UTV) and that there 
is a Haar-distributed random unitary matrix V, independent of U, such that for each j, 
Mj = VMjV*, with Ml,..., a collection of deterministic matrices also satisfying (2). 


3.2. Tails of the series. Let us hrst prove that Theorem 2.4 can be deduced from the 
particular case where there is hq such that for all n, we have 

\n\ > no Vj = 1,..., fc, ttnifj) = 0. 

Let £ G (0, a/2) such that the domain of each fj contains the annulus of complex numbers 
z such that a — 2e < \z\ < b + 2e. 

Lemma 3.1. There is a constant C independent of N such that for any n such that < N 
and any j = 1,..., /c, we have 

E I Tr A”M,f < Cn^ (l„>o(6 + + ln<o(a - 


Proof. With the notation of Section 3.1, let E v denote the expectation with respect to the 
randomness of V. For each n G Z and each j, by Lemma A.6, we have 

Ev|TrA”Mjf = EvTr A^VM^V* Tr(A*)”VM*V* 


A2 


1 


Tr A^ Tr(A*)^ Tr Tr M* + Tr A^(A*)" Tr 


Tr A'^ Tr(A*)’" Tr M^M* + Tr A* (A*)” Tr Tr M( 


< 


A(A2 - 1) 

1 Tr A"|2 Tr Tr M* + Tr A"(A*)’" Tr 


A2 - 1 

< C(|TrA’*p +A-^TrA^(A*)’^) , 

where C* is a constant independent of N. Then the conclusion follows from Lemma A. 5. □ 

Lemma 3.2. There are some constants C > 0 and c G (0,1) and a sequence E = En of 
events such that 

P(T) ^ 1 

N^oo 

and such that for all N, all ni > 0 and all j = 1,. .., k, we have 

E a„(/,)Tr(A-M,)| < ^^(l - c)"L 

|n|>ni 


Proof. By [16, Lem. 3.2], we known that there is a constant Ci such that the event 
E = En-.= {Vn > 0, II A^llop <Ci{h + £)"} n {Vn < 0, || A"||op < C'i(a - e)^} 
has probability tending to one. 







12 


FLORENT BENAYCH-GEORGES AND JEAN ROCHET 


Then one concludes easily, noting first that by non-commutative Holder inequalities (see 
[2, Eq. (A. 13)]), we have 


{ C'i(6 + e)'^N^N-^ Tr M^M* if n > 0 

Cl (a - e)^N^N-^ TrMjM* if n < 0 

and secondly that there is c G (0,1) such that for each j, the sequences 


^nifj ) 


{b + ey 

(W 


n>0 


^nifj ) 


(a — ey 

(T^ 


n<0 


are bounded. 


□ 


As a consequence of Lemmas 3.1 and 3.2, for any 0 < no < ni < and any j = 1,... ,k, 


E 


If E «»(/j)Tt(A”M,) 

|n|>no 


< 


no<|n|<ni 

+e|i£ a„(/,)Tr(A”M, 

|n|>ni 


< C\an{fj)\n‘^ {in>o{b + sY"' + ln<o{a — ey"''j 

no<|n|<ni 

+CN{1 - 

Choosing first ni = [AlogAj for A a large enough constant and then using the fact that 
for any j = 1,..., k, 


^ |a„(/j)|n^ {ln>o{b + + ln<o(a - < oo, 

nGZ 


we deduce that for any 5 > 0, there is no > 0 fixed such that for all N large enough. 


5^E a,,(/,)Tr(A-M,) 

i=l |n|>no 


< 5, 


( 7 ) 


for jv as in Lemma 3.2. 


Let us now suppose Theorem 2.4 to be proved in the particular case where there is no such 
that for all n, we have 

|n| > no Vj = 1,..., fc, ttnifj) = 0 

and let us prove it in the general case. Let X^r denote the random vector of (3). We want to 
prove that as A —)■ oo, the distribution of tends to the one of Q := ■ ■ ■, G{fk)), 

i.e. that for any function F : —)■ C which is 1 Lipschitz and bounded by 1, we have 
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J=1 


To do so, we first set 

XN,no ■■= ( a„(/,)Tr(A-M,)-Tr(M,)ao(/, 

|n|<no 

By hypothesis, for any hxed rio, converges in distribntion to a centered complex 

Gaussian vector g^o ■= {Qnoifi), • • •, Qnoifk)) such that 

n<3nAh)<3n,{f,)] = n<3Ui)<3U,)] + t 
= E 


-no 

lij 


E 




6(/i)<?(/i) +C. 


where lim lhl^° I + 1*^11° I = 0- Therefore, 

riQ^cxD ^ ^ 

1<2 j</i; 

|E|F(A'A,)-F(e)|| 

< |E |F(A'„) - F(A«^„.)|| + |E 1F(Aa,,„.) - F(ftJ]| + |E |F(6„.) - F(a)|| 

< 2P(F') + E [ijII-Ya, - Ak,„ji^] + |E|F(A>„,) - F(e„„)ii + E [iie„, - aiig 

which can be as small as we want by (7) and the fact that X 7 v,no -^Qno if Gno and Q are 
coupled the right way. 


3.3. Proof of Theorem 2.4 when the /j’s are polynomial in z and z~^. We suppose 
here that there is no > 0 such that for all n > tiq and all I < j < k, an{fj) = 0. Without 
any loss of generality, we also assume that for all j, ao(/j) = 0. In this case, any linear 
combination of the Tr/j(A)Mj’s can be written 

k 

Gn := 5^n,Tr/,(A)M, = TrA"N„ 

i=l \n\<nQ 


k 

where N„ := z/ja„(/j)Mj. Written this way, we notice that to prove that the limit 
i=i 

distribution of Gat is Gaussian, we simply have to prove that the random vector 


converges in distribution to a Gaussian vector. We will prove it by computing the limit of 
the joint moments. 


Before going any further, recall that by the preliminary randomization of the N^’s from 
section 3.1, we suppose that A = UT (instead of A = UTV) and that there is a Haar- 
distributed random unitary matrix V, independent of U, such that for each j, Nj = 
VNjV*, with Nj a deterministic matrix. 

We shall proceed in three steps: 
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a) First, we prove the asymptotic independence of the random vectors 

(TrA"N„,TrA-"N_„), 


' n>l 


by the factorization of the joint moments. More precisely, we prove, thanks to 
Corollary A.4, that for any (pn)nLi, (gn)”ii, (sn)nLi, 


E ^ Y[ ( Tr A’^N^) ^" ( Tr A«N„) (Tr A-"N_„) "■" (Tr A-"N_„)' 

l<n<no 


= JJ E (Tr A”N„)^"(TrA’^N„)''"(TrA-"N_„)’'"(TrA-^N_„)'" + o(l) 

l<n<no 

b) Then, we prove for any hxed n, the random complex vector 

(TrA’*N„,TrA-"N_„) 


converges in distribution to a centered complex Gaussian vector thanks to the 
criterion provided by the Lemma A.7. This criterion consists in proving that the 
joint moments, at the large N limit, satisfy the same induction relation as the 
moments of a Gaussian distribution. 

c) It will follow from a) and b) that when all f/s are polynomials in z and z~^, the 
random vector of (3) converges in distribution to a centered Gaussian vector. To 
conclude the proof, the last step will be to prove that the limit covariance is the 
one given in Theorem 2.4. 


In the proofs of a) and b), we shall need to compute expectations with respect to the 
randomness of the Haar-distributed matrix U. More precisely, we shall need to compute 
sums of expectations with respect of U resulting from the expansion of products of traces 
involving powers of A (such as Tr A"’N„). To do so, we will use the Weingarten calculus 
(see Proposition A.l) and shall always proceed in the following way: hrst, we use (28) to 
state that all the terms of the sum are null except those for which the left (resp. right) 
indices involved in u are obtained by permuting the left (resp. right) ones involved in u. 
Then, we claim, by Remark A.2, that among the remaining terms, we can neglect all those 
whose indices are not pairwise distinct. At last, once all the remaining terms are, up to 
multiplicative constant, equal to Wg((T) for some permutation a, we neglect all those for 
which a ^ id (see Remark A.3) and the summation hnally gets easy to compute. We 
introduce here a notation that we shall use several times : 

If := {(!i,.... 4 )e{i....,ivr; R, ... ,in are pairwise distinct} (8) 

(this set implicitly depends on N). 


3.4. Proof of b). In this part, as n is hxed, we shall denote N„ (resp. N_„) by M = [M^] 
(resp. K = [Kij]). For any non-negative integers p, q, r, s, wet set 

mp,q,r,s ■= E ( Tr A’^M)^(Tr A^M)"( Tr A-’^K) ( Tr A-^K)' 
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and our goal is to show that, as N goes to inhnity, the numbers mp^g^r,s have limits satisfying 
conditions (34), (35), (37), (38) and (39) of the Lemma A.7. Note that (34) and (37) follow 
from the fact the the Haar measure on the unitary group is invariant by multiplication by 
any e^®, 0 G M. We shall use the following notations 


lim — Tr (MM*) =: Um ; 
N^OO N 


lim — Tr M =: tm 

N^oo N 


lim — Tr (KK*) =: Uk 

lim — Tr K =; tk- 

N^oo N 


lim ^ Tr (MK) =: /3 m,k 

N^OO TV ^ 


3.4.1. Tr A’^M and Tr A ”'K are asymptotically two circular Gaussian complex variables 
satisfying conditions (35) and (39). We simply have to show that for any integer p > 1 

E|TrA"M|'^ = p\{b^^{{n-l)\TM\^ + aM)Y + o{l) , (9) 

E|TrA-"K|"'’ = p!(a-2"((n-l)|rK|' + aK))'’ + o(l) , (10) 

E(TrA"M TrA-”K)^ = p! ((n - l)rMrK +/3m,k)'’+ o(l) . (11) 


We shall prove it by induction on p. So hrst, we show the previous relation for p = 1. Recall 
that we assume that M = VMV* and K = VKV*, where M and K are deterministic, so 
that, using the Lemma A.6, we have (denoting again by Ev the expectation with respect 
to the randomness of V), 


Ev TrA'^VMV* 


Ev TrA-"VKV* 


EvTrA^VMV TrA'^VKY* 


^TrA’"(A*)” (^^TrMM* - |^TrM|^^ 

iivA”riiivMr+o(i) 

f IV A-"(A-)-" (i IV (KK-) - I ^ TV' K 

I TrA-’^l^l —TrKlVof—^ 

I I liv I \NJ 

— TrMK - — TrM—TrK 
N N N 

TrA”TrA-”^TrM^TrK + o|^^^ . 


This is asymptotically determined by the limits of E | Tr A*^ | ^, E | Tr A ” | ^, E Tr A"' Tr A , 

A“^E Tr A"'(A*)"' and A“^E Tr A“"(A*)“"'. First, we compute E | Tr A""!^ for n > 1. We 
write 

E I Tr A I = ^ ) IE [uijjj ‘ ‘' ‘ ‘' '^jnji] A’l ''' (12) 
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From (28), we have a condition on the i^s and the j^s for a non-vanishing expectation, 
which is the multiset^ equality 

= {jli ■ ■ ■ ijn}mi (13) 

The first consequence of (13) is that the sum is in fact over 0{N^) terms which all are at 
most 0{N~^), which means that any sub-summation over o(iV") terms might be neglected. 
So from now on, we shall only sum over the n-tuples (ii,..., i„) G (recall notation (8)). 
Then (12) becomes 
|2 


E Tr A" 


E 


2 2 

Sir • • Si, 


E® 

CT&Sn 


‘ ' ' (n)*a(l) 


+ 0(1) 


Let c G S'„ be the cycle (12 ■ ■ - n). From (28) (see Remark A.2), as long as the 4+ are 
pairwise distinct, one can write 


E 




= Wg (cTC ^cr ^c) 


and from (29) and Remark A.3, we know that the non-negligible terms are the ones such 
that ac~^(j~^c = id, i.e. ac = ca, which means that a must be a power of c and so, 
only n permutations cr contribute to the non negligible terms. At last, as Moeb(i(i) = 1, 
we have 

r = X nA-”(l + o(l))+o(l) 

(il 5 • ■ ■ Jn ) 

N 


E Tr A" 


= n \ — 
N 


E 

2=1 


s" I -h o(l) = -h o(l). 


In the same way, one can get 

E|TrA-”|^ = na-2’^ + o(l), 
ETrA'^TrA"^ = n + o(l). 

Let us now consider A“^E Tr A"'(A*)"' for n > 1. We have 
1 


—ETrA'^(A*)" = A-i^E[ 


Jiioil ■ ■ ■ 


' u. 


1—ijn] 


SjiSji 


' Si^S 




(14) 


l<20,2l,...Jn<Y 
’^0—Jo? in—jn 


As previously, we know that by (28), that for the expectation to be non zero, we must 
have the multiset equality 

{* 0 , ■ ■ ■ ~ {jo-} ■ ■ ■ ijn}mi ( 13 ) 

The first consequence of (15) is that the sum is in fact over terms which are all at 

most 0{N~^~^), so that any sub-summation over terms might be neglected. As 


^We use the index to in { • }^ to denote a multiset, which means that {xi,..., Xn}^ is the class of the 
n-tuple {xi ,..., Xn) under the action of the symmetric group Sn- 
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previously, we shall sum over the pairwise distinct indices (see notation (8)). Hence 
(14) becomes 


AT-ie TrA’"(A* 


N-' 


2 2 

^*1 ■ ■ • 


(T^Sn+1 

o-(0)=0 

cr(n)=n 




' U. 


^cr(n —l)^cr(n) 


Let c G S'„+i be the cycle (0 12 ■ ■ -n). From (28) (see Remark A.2) one can write 


E 


TIzqU ^^71 —l^n^G(0)L(l) ^G(n —l)G(n) 


= Wg( 


ac V 


)• 


As previously, a must be a power of c for the asymptotic contribution to be non-negligible. 
However, this time, we impose (t( 0) = 0 and <j{n) = n, so that the only possible choice is 
a = id which means that only one term contributes this time. At last. 


—ETrA’"(A*)” = b^^ + o{l), 


The same way, one can get 


^E TrA-’"(A*)-^ = a-^^ + o{l). 

This concludes the first step of the induction. 

In the second step, we have to prove the following induction relation: for any p> 2, 

E|TrA"M|^^ = pE| TrA”M|^E| Tr + o(l) (16) 

EjTrA-’^Kl^^ = pE| TrA-”K|^E| Tr +o(l) (17) 

E (Tr A"M Tr A-”K)P = pE [Tr A"M Tr A'^K] E [(Tr A"M Tr A'^K)^"^] 

+0(1) (18) 

I 1 

Let us hrst consider E | Tr A^'M] . We shall use the following notation 

Tr A M ^ ^ ^ ^ , 

where the bold letter i denotes the (n + l)-tuple (zq, ■ ■ ■ ,in) and where we set 


•- ^20^1 * * ‘ ^2rj —f2 


Si-^ • • • Si^ . 


(19) 


Hence, 


E I Tr A’"M| ^ = ^E [uii • • ■ UipUji • • • mjp] ■ ■ ■ SipMiPiPSjpMjPjP (20) 


ji iP 

1 ,.,1 


As usual, we know we can sum over the i^’s satisfying that (i^,..., F) (the p{n + l)-tuple 
obtained by concatenation of the i’s) has pairwise distinct entries and such that we have 
the set equality: 

{ij, l<A<p. 0<;i<n} = {j^, l<A<p, 0<p<ii}. 


( 21 ) 
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Then, in order to have Wg{id), we must match each of the {n + l)-tuples ..., with 
one of the (n + l)-tuples ..., i.e. that for all 1 < A < p, there is a 1 < A' < p such 
that we have the set equality 



Iq, ll, 



JO ’)j 1 ? 




We rewrite (20) by summing according the possible choice to match {i^} = ... Jh} 


E I Tr A"M| ^ ^ 5^ E [uii ■ ■ • UipUji • • • up] • • • s^MiV^ivs-^pMjVjV + o(l), 

^=1 


where stands for the set equality {zq, i\, ■ ■ ■, ih} = ■ ■ ■ ,jn}- Then, we know 

that the set of indices {zg, • • •, ih} is disjoint from the others, so that by Corollary A.4, 

E [uii ■ ■ ■ UipUji ■ ■ ■ up] = E [uiiupT] E ^Ui2 ■ ■ ■ UipUp ■ ■ ■ up-iup+i ■ ■ ■ up] 


and up to a proper relabeling of the indices, all the choices lead to the same value of the 
expectation, so that 


E Tr A’^M 


2p 


p ^E[ni2' 

(j V 


• Mipl S\2Mj2flS\2M^2j2 ■ ■ ■ S\P MjP jP S\P MaP ^P + o(l) 
J J 1 J JnJO ’‘nf'Q J JnjQ V / 


* *=^n+l 

J ^-^n + l 


pE 


Tr(A’"M)| 


E 


Tr(A”M) 


2(p-l)' 


+ 0 ( 1 ). 


This proves (16). In the same way, we prove (17) and (18), and thus conclude the proof of 
the induction. 


Remark 3.3. In the last computation, we split the expectation and so we separated the 
summation implying that 


p(n+l) 


— I^ X I^ 

- ++1 -L(p_i)(n+1) 


which is obviously inaccurate. Nevertheless, we easily see that 


CardlJ^^,) = Card (^I^+i x (l + o(l)), 

which means that this inaccuracy is actually contained in the o(l). 


To conclude the proof of b), we have to prove that Tr A"M and Tr A ""K satisfy Condition 
(38) at the large N limit. 












FLUCTUATIONS FOR ANALYTIC TEST FUNCTIONS IN THE SINGLE RING THEOREM 


19 


3.4.2. Tr A"'M and Tr A“”K satisfy Condition (38) at the large N limit. We apply the 
same idea as previously, but for a slightly more complicated expectation. Let p, g, r, s be 
positive integers and such that p — q = r — s. We denote joint moments by mp^q^r,s- 


and as 


m 


p,q,r,s 


TrA”M 

TrA-’^K 


E(TrA"M)^(TrA^M) ( Tr A-”K)'’( Tr A-’^K) , (22) 




^ ^ '^in-ii-n^i„ ' ' ''^ioh^h^io'in ~ ^ ^ ^ioinl 


lv!^0 


we rewrite (22) as follow 


E E If 



1<A<P 

l</i<g 

l<u<r 

i<e<s 






ALa MKi^v i^i^K£ e £B UixUieU^fiU]^!' 




(23) 


(recall that the Si = Sj^ • • • for i = (io, • • •, in))- As previously, we deduce from Propo¬ 
sition A.l that for the non vanishing expectations, we must have the following multiset 
equality 

{ip, l<A<p,0</i< n}^ U l<A<s,0<p< n]^ 

= {j'i, l<A<g,0</r< n]^ [j [kp, l<A<r,0</r< n}^ , (24) 

from which we deduce that we can restrict the summation to the tuples such that 




and that, for the non negligible terms (he. those which lead to Wg{id)), we must match 
each of the {n + l)-tuples involved in u (the i’s and the I’s) with one of those involved in 
u (the j’s and the k’s). For example, we sum according to the choice the “partner” of i^. 


m. 


•p,q^r,s 




“=i (i,i)eP i<A<p 
(j.k)eP llOl? 




1 l<iy<r 

i l<0<s 


ee E n 

/3=i (i,i)eP 1<^<P 
(j,k)eP 


5jA5j/j' ^ ^ _ _ 




^e-UiA-Up-UjA.-Uk'' + o(l) 


Jink'S 

1<0<S 


where, to simplify, (i, 1 ) stands for the (p + s)(n + l)-tuple obtained by the concatenation 
of the i'^’s and the F’s, and implicitly means As previously, we use the 
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Corollary A.4 to split the expectations. Hence, one easily gets 


m. 


'p^q,r,s 


= qE 


Tr A”M 


mp_i^q_i^r,s + rE [Tr A”M Tr A ’^K] mp_i^q^r-i,s + o(l) 


To get the other relations, we just sum according to the choice of the partner of (resp. 
and T). 


3.5. Proof of a): asymptotic factorisation of joint moments. The proof relies mostly 
on Corollary A.4. We hrst expand the expectation 


E [ Yl ( Tr A’^N^) ^" (Tr A”N„) (Tr A-”N_„) "■" ( Tr A-"N_„) 

l<n<no 


(25) 


Let denote the (i, j)-th entry of N„ and recall that for i = (io, ■ ■ ■ ,in), we set 


tl'i . UiQii ' ' ' , 5; . 5^^ ' ' ' Si 


We get 


TrA"N„ = 






(n) 

^0 




— 1 ? • • • 


so that 


^ H Yl n 


l<n<no l<\<pn 

jn,l jn,gn 
T-n,l \^n,rn 


c c ^ Jo in' C 

S\^n,u S\n,e ^0 ^ 0 n 


where we use bold letters such as i”’'^ to denote (n + l)-tuples (zq’'^, ... ,z”’^). We can 

use the same ideas as in [16, Lemma 5.8] to state that the non-negligible terms of the sum 
must satisfy that for all n, there are as much [n + l)-tuples involved in u as in u, which 
means that 

Pn T Qn T i’ri) 

and that we must have the multiset equalities 

no 

1< A<p„,0</z<n}^U{C^ 1< A<s„,0</i<n}^ 

n=l 

no 

= \J{r/, l<X<qn,0<p<n}^[j{k;’\ l<A<r„,0</x<n}^. (27) 

n=l 

We deduce that there are a non zero terms in (26) and we can easily show 

that any subsum over a is negligible so that for now on we shall sum over 
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the non pairwise indices. Then, we know that we can neglect any expectation E u which 
won’t lead to Wg{id) (see (29)) so that (27) becomes 

VI < n < no, 1 < A < 0 < /I < n} U 1 < A < 0 < < n} 

= 1 < A < gn, 0 < /i < n} U 1 < A < r„, 0 < /i < n} . 

It follows that the set of indices involved in the expansion of the Tr A"'N„, TrA“"'N_„, 
TrA’^N^, TrA-’^N.^, is disjoint from the set of indices involved in the expansion of the 
Tr A"^Nm, Tr A""*N_m, Tr A^^N^, Tr A-"*N_m, as long as n 7 ^ m. Therefore, Corollary 
A.4 allows to conclude the proof of a). 


3.6. Proof of c): computation of the limit covariance. Let f,g be polynomials in z 
and z~^ and let M, N be A x A deterministic matrices such that, as A —>■ oo, 

^TrM ^ T ; ^TrN ^ r ; ^TrMN* ^ a ; ^TrMN ^ /3. 

A A A A 

We need to check that the limits of both sequences 

E (Tr/(A)M - aoif) TrM)(Trg(A)N - ao(g) Tr N) 

and 

E (Tr/(A)M - aoif) Tr M)(Trg(A)N - aoig) Tr N) 

are the ones given in the statement of Theorem 2.4. Note that it suffices to compute the 
limits for f = g and M = N. Indeed, using the classical polarization identities for M and 
N, first for general polynomials /, g, we reduce the problem to the case M = N. Then, we 
use polarization identities again to reduce the problem to f = g. 


Also, recall that since A A for any deterministic 6, we know that for any positive 
distinct integers p, q, we have 


E TrAPMTrA-«M = E Tr A^’MTr A-iM = 0. 

It follows, using (9), (10) and (11), that 

E |Tr/(A)M-ao(/)TrM|2 = a^(/)M7)E TrA™MTrA-M 

7^0 


E(Tr/(A)M-ao(/)TrM)2 = 


(|a„(/)|"E I Tr A”M|' + |a_„(/)|2E | Tr A-"M|") 

n>l 

+ |a_„(/)|^a"^”) ((n - l)|r|^ + a), 

n>l 

a^(/)an(/)ETrA™MTrA"M 

7^0 


= 5^2a„(/)a_„(/)E Tr A^MTr A'^M 

n>l 
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which concludes the proof. 


2 «n(/)a-n(/) ((^ - 1)1-^ + , 

n>l 


4. Proof of Corollary 2.14 


with Af := ^ if |z| < a. 

withSf :=-TrE.>iSif|^l>^ 

(in the first case, we used the fact that | det A| = det T). Then, by Theorem 2.4, 

converges, for the finite-dimensional distributions, to a centered complex Gaussian process 

i'^z)\z\<a U {^z)\z\>b 

with covariance defined by 

_ 

EA^Az' = 0, EAzAz' = -log(l- -), 

EBzBz' = 0, EB^Bz'=-log{l-^=), 

zz' 

EAzBz'=-\og{l-—), EAzBz' = ^, 

z 

where log denotes the canonical complex log on i?(l, 1). Then, one concludes by noting 
that for A^B G C, 2 Re A Re i? = Re(Ril -|- AB). 

Appendix A. 

A.l. Weingarten calculus on the unitary group. Here we summarize the results we 
need about integration with respect to the Haar measure on unitary group, (see [25, Cor. 
2.4 and Cor. 2.7]). 

Let Moeb denote the Mobius function of the lattice of non-crossing partition, defined for 
example in [44, Lect. 10]. To simplify, for any fc-tuples i = (R,..., 4) and j = (ji,... ,jk), 
we set 

“ij := 


It is easy to see that, for any z ^ S, we have 
log I det( 2 ; — A)| = 


TrlogT + ReAf, 

;7V 
z 


N\og\z\+ReB^, 


'^iljAi2j2 ■ ■ ■ '^ikjk 
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Proposition A.l. Let k be a positive integer and U = {uij) a N x N Haar-distributed 
matrix. Let i = (zi, ... ,4), i' = (4,... ,4), j = (ji,... ,4) and j' = {j [,..., be four 
k-tuples of {1,..., N}. Then 


a,T£Sk 


( 28 ) 


where Wg is a function called the Weingarten function. Moreover, for a G Sk, the asymp¬ 
totical behavior ofWg{a) is given by 

iV'=+Hwg(a) = Moeb(a) + o(^^^, (29) 

where |(t| denotes the minimal number of factors necessary to write a as a product of 
transpositions. 


Remark A.2. One should notice that if all fc-tuples (4, • • • ,ik), (ji, • • • ,jk), {i'l, ■ ■ ■ ,i'k)j 
and {j[,... have pairwise distinct entries, then (28) becomes simpler because in this 
case there is at most one non-zero term in the sum. 


Remark A. 3. The permutation a for which Wg((j) will have the largest order is the only 
one satisfying |cr| = 0, i.e. a = id. Also, Moeb(z(i) = 1 (see [25]). 


Here is a useful corollary which permits to simplify many computations. 

Corollary A.4. Let i = (4,..., 4 ), j = (ji,...,jp), k = {ki,...,kg), 1 = (4,..., 4 ), 
i' = (4,.. ■, ip), y = ifi, ■ ■ -Jp), k' = {k[, ..., kg), V = (4, ...,i'g) be tuples such that the 
multisets defined by i and i' (resp. by j and}', by k and k', by 1 and V) are egual and such 
that 

{ii,...,ip}n{ki,...,kg} = {ji,...,jp}n{li,...,lg} = 0. 

Then 

E [uijMk.iMi'j'MkM'] = E [uijWj] E [Mk,i«^] (^1 + ^ (7^) ) • 


Proof. To prove this result, we hrst recall the exact expression of the Mobius function : for 
any permutation a with cycle decomposition C'iC '2 ■ ■ ■ Cr, 

r 

Moeb(cr) = ]^(-l)l'^‘l“^ Cat|Ci|-i, 

i=l 


where Cat^ is the fc-th Catalan number. 


2k 


k 1\ k 

permutations with disjoint supports, then 

Moeb(cr o r) = Moeb(cr) Moeb(r) and 

Thus 

jgp+q+\(TOT\ \Yg(cr or) = Moeb(cr o r) -\- O ^ 


. Then, obviously, if cr and r are two 


cr o r = cr N- r . 
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Moeb((T) Moeb(r) + O ( j 


= n 


P+l-l Wg(a)Ar'?+l"l Wg(r) (^1 + 0 


So that 

Wg(cr o r) 

One can easily conclude. 


Wg(a)Wg(r)(l + 0( — 


□ 


We also need the following lemmas in the paper. 

Lemma A.5. Let A = UT with U Haar-distributed on the unitary group and T deter¬ 
ministic. Let £ > 0. There is a finite constant C depending only on e (in particular, 
independent of N and ofT) such that for all positive integer n such that n® < (2 — e)N, 
we have 

( 2 \ ^ 

m2 + j (30) 

and 

E[|TrA-|2] < c(m2 + '^'^ , (31) 

where m 2 := A“^TrTT* and m^o ■= ||T||op. 

Proof See [11, Th. 1]. □ 

Lemma A.6. Let V be Haar-distributed and let A, B, C, D be deterministic N x N ma¬ 
trices. Then we have 

E TrAVBV*TrCVDV* = —^—(Tr A Tr C Tr B Tr D + Tr AC Tr BD) 

— 1 

“A(A^—^(TrATrCTrBD + TrACTrBTrD) 


Proof. Let AlAr(C) denote the set of A x A complex matrices. It has already been proved, 
in [16, Lem. 5.9], that for any matrices A, B, C, D G AlAr(C), we have 


E Tr AVBV*CVDV* 


^ {TrACTrBTrD + TrATrCTrBD} (32) 
- —- {Tr A C Tr B D + Tr A Tr C Tr B Tr D } . 


We deduce that 
EV(g)V* 0 V® V* 


1 

A2 - 1 



E. 




Ep,i + E, 




Ek,e ® Ei^k ® E^ 
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~ ^ ^ Ej,k ^ Efj 0 Ei^£ ® Ek^i + Eij ® Ej^k ^ Ek^t ® Ef^i j , (33) 

where the Er^s denote the elementary matrices. Indeed, the linear morphism \h from 
A^ 7 v(C)®^ to the space of 4-linear forms on A4n{C) defined by 

^(M^N(g)P^Q)(A,B,C,D) := TrAMBNCPDQ 

is an isomorphism, and (32) proves that the left and right hand terms of (33) have the 
same image by 'h. Then, applying 

M(g)N(g)P(g)Qi— >Tr (AMBN) Tr (CPDQ) 

on both sides of (33), we deduce the lemma. □ 


A.2. Moments of a Gaussian vector with values in C^. The following lemma allows 
to prove that a complex random vector {X, Y) is Gaussian without having to compute all 
its joint moments, by only proving an induction relation. 


Lemma A. 7. Let {X, Y) be a Gaussian random vector with values in whose distribution 
is characterized by 


EX 

= EX^ = EX = EX^ = 0 , 

(34) 

and 



E|Xp = ax ; 

E|X |2 = ay ; EXX = axv ; EXX = 0 . 

(35) 

Then, the moments 

■■= EX^’XV^X" 

(36) 

satisfy 



P 

— q ^ r — s fnp^q^r,s = 0, 

(37) 


^ 1,5—1,r,s “1“ l,5,r—l,s 

_ 1 1,5—l,r,s “1“ ^^Xy^p,5—l,r,s—1 

1 5f^y^p,5,r—l,s—1 “1“ P^Xy^p—l,5,r—l,s 

^^y^p,5,r_l,s_l “1“ ?^yy^p,5—l,r,s—1 

(38) 

and 

mp,o,p,o = E(XX)^ = p!a^y. 

(39) 


Conversely, if (X, Y) is a random vector with values in such that both X and Y are 
Gaussian and have joint moments rhp^q^r,s satisfying (34), (35), (37), (38) and (39), then 
(X, Y) is Gaussian. 





26 


FLORENT BENAYCH-GEORGES AND JEAN ROCHET 


Proof. First, one easily obtains (37) by noticing that for any 6* G M, (e'^X, e (X, Y). 

To prove the remaining, we consider (X, Y) as a real 4-tuple (3fJ(X), 7y(X), 3fJ(X), 7J(X)) =: 
{xi,X2,X3,Xi) with covariance matrix 

/ax 0 ^(axy) ^{(^xy) \ 

0 ax ^(axy) -^{(^xy) 

^{axy) ^{axy) ay 0 

\^^5(crxy) -^{axy) 0 ay J 

Its Fourier transform is given, for t = (ti, ^ 2 , ^ 3 ,^ 4 ), by 

$(t) := E exp (i(fia:i -1- 12X2 + hx^ + 14X4)) 

= exp { - ^{ax{tl + tl) + ay{tl + tl)) - ^3fJ(axy)(ti^s - ht^) - ^^(crxy)(^ 2^3 + 


F := - 
2 


We dehne the differential operators 

+ 1^2 

dy = 83 + 184 


8x = 81- 182 

8 y = 83 — 184 


so that 


and 


Ex^’xVw' = 


'9x4’(t) — — -(crx(ti + 1 ^ 2 ) + c’‘xy(^3 ~ i^4))*f’(t) 

= — 2 — 1^2) + <Txy(t3 + h 4 ))‘h(t) 

9y<h(t) = -^(cry(t3 + 1 ^ 4 ) + <Txy(tl - h 2 ))$(t) 

9y$(t) = -]-[ay{t3-it4)+ aYV{ti+it2))^{t) 


(40) 


We easily deduce that 

- roxyS;“'9pyaf4(t)|,^„ 

so that 

EXPXVw' = gcrxEX^’-iX''“V’’F'+ rcrxyEX^’-iX‘'W-iF'. 

This proves (38). To prove the last point, we simply write 

E{XYY = {-iy8Px8^^{t)l^^ ; «F(t) = ( - ^ +^^xy{t3 - 
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Then, using Leibniz formula and noticing that for all k < p, 

Qk + 1^2) + ~ 1^4) 

one can easily conclude. 

Conversely, let {X, Y) be a random vector with values in such that both X and Y are 
Gaussian and have joint moments rhp^q^r,s satisfying (34), (35), (37), (38) and (39). Let N 
denote the set of non-negative integers and let us define the sets 

/C := {(p,g,r,s) e p - q ^ r - s] , 

'kLo := |(p, g,r, s) e ; {r + s + \p-q\) {p + q+\r - s\) {p + r+ \q-s\) {q +s + \p-r\) = 0 
and, for A; > 1, 

kik ■= {(p, g,r, s) e (p - 1, g - 1, r, s) or (p - 1, g, r - 1, s) or 

{p,q- l,r,s - 1) or ip,q,r - l,s - 1) G kik-i} 

Then by hypothesis, the joint moments function fhp^q^r,s coincides with mp^q^r,s on /C U "Ho- 
Besides, by (38), if fhp^q^r,s coincides with mp^q^r,s on kik-i, then fhp^q^r,s coincides with 
mp^q,r,s on 'Hfc. As 

N" = /CuljTffc, 

k>0 

we deduce that fhp^q^r,s coincides with mp^q^r,s on wich implies that {X,Y) (X, Y). □ 

Acknowledgments. The authors wish to thank the anonymous referee for his careful 
reading and his useful advices. 
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